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Abstract 

IGMtransmissiorj^ is a Java graphical user interface that implements Monte Carlo sim- 
ulations to compute the corrections to colours of high-redshift galaxies due to intergalactic 
attenuation based on current models of the Intergalactic Medium. The effects of absorption 
due to neutral hydrogen are considered, with particular attention to the stochastic effects of 
Lyman Limit Systems. 
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1 Introduction 



IGMtrcinsmission is a Java graphical user interface (GUI) that models the effect of intergalactic 
neutral hydrogen attenuation on the colours of objects based on the model of Meiksin(2006). The 
code allows the effects to be modelled for a range of galaxies and filters provided. The galaxy spectra 
are taken from the stellar + nebular emission line models of Leitherer et al. (1999). Photometric 
filters are included for the Hubble Space Telescop^ the Keck telescop^ the Mt. Palomar 2OO-incl]0 
the SUBARU telescop^ and UKIDS^ Alternative spectra and filters may be straightforwardly 
added by the user. 

To model the colours of high-redshift galaxies a model galaxy spectrum is taken, shifted appro- 
priately for a chosen redshift, and then attenuated to account for intervening neutral hydrogen using 
then chosen model and set of parameters. Monte Carlo simulations over many lines of sight are 
performed, with the option of obtaining each colour independently or calculating an average. To do 
this, a population of discrete, Lyman Limit Systems are drawn from the distributions, dN/dz and 
dN/dri, outlined in Section 2.1, placed along random lines of sight and their total contribution to 
the intergalactic transmission function calculated. In addition to the default choice, the alternative 
Lyman Limit Systems distribution of Inoue & Iwata (2008) is allowed for. Alternatively, the pa- 
rameter values defining the model may be chosen by the user. Specifying a fixed set of individual 
Lyman Limit Systems is also accommodated. Mean contributions to the opacity due the Lya forest 
and the photoelectric absorption of the optically thin IGM are then included. The colours and IGM 
k-corrections are computed through the chosen broad-band filters, given this amount of attenuation 
for many lines of sight. The cumulative probability distributions, over the different lines of sight, are 
also produced. The model spectra, broad-band filter combinations, attenuation model parameters 
and output options are all chosen via the GUI. 

The code requires Java v. 1.6. It has been tested on Windows Vista'^^Home Basic, Debian 
Linux v5.0 and OS X v. 10.5.8. The code uses the random number generator included in the Sun 
Microsystems' Java Development Kit v. 1.6. To query the default version of Java running on a given 
system, type "Java -version." If it is not v. 1.6, assistance from a system administrator may be 
required to obtain it. A screenshot of the GUI and an explanation of how it is used is included in 
Appendix |Xj The details of the code design are provided in the next section, followed by a section 
describing the various code features. 



2 Code design 

2.1 Modelling intergalactic attenuation 

Intervening absorption systems are likely to be of various origins; for example the outer edge of 
galaxies, halo gas and diffuse material in intergalactic space. The absorbers are generally considered 
to be discrete and separate systems. They are classified into three main categories depending on 
their properties: Lya forest absorbers (with column densities Nm ^ lO^^cm"^), Lyman Limit Sys- 
tems (LLSs) (lOi^cm-2 < Nm ^ lO^^cm-^) and Damped Lya Absorbers (DLAs) (A^'m ~ lO^^cm-^) 
(Meiksin 2009). It has been suggested that there is an upper-cut off at 3 — 5 x 10^^ cm~^, and systems 
with column densities lower than lO^^cm"^ may exist but are difficult to detect. The Lya forest 
absorbers are by far the most common, and the DLAs the rarest and most easily identifiable in spec- 
tra. The intermediate LLSs have a sufficiently large column density to absorb photons with energies 

■^http: / /hubblesite.org 

^http:/ /www. keckobservatory.org 

^http: / / www. astro. caltech.edu/palomar 

^http: / / www.naoj .org 

'^http; / / www.jwst.nasa.gov 
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above the photoelectric edge (the Lyman hmit). The classifications are not strictly exclusive; for 
example, DLAs will also produce LLSs. The total effective optical depth of the IGM is the sum of the 
optical depths due to each absorbing system. The two main effects responsible for attenuation of the 
spectrum of a distant light source are the photoelectric absorption by LLSs and optically thin IGM 
and resonance line scattering by the Lja forest. Additional contributions are made by intervening 
metal systems and helium. Contributions of metals and Hel (neutral helium) are small (Madau 
1995), while Hell (singly ionised hehum) will contribute only at wavelengths A < 228(1 + z) A. This 
is outside the wavelength range considered. Due to the discrete and clumpy nature of these systems, 
it is not possible to know the amount of absorption along one particular line of sight. However, it 
is possible to predict a mean amount of absorption and model the statistical fluctuations. Madau 
(1995) uses analytical approximations to derive the mean attenuation. More recently work has been 
done using Monte Carlo simulations to predict the mean attenuation for the models used in the code. 



Lyman transitions 

The Lya forest does not produce significant fiuctuations, therefore a mean contribution is assumed. 
The mean optical depth for the Lyman transition n — > 1 is defined by 

fn = -ln(exp(-r„)) (2.1) 

where (exp(— r„)) is the corresponding mean transmitted flux and n corresponds to a, /5 etc. The 
values of r„ used in this investigation for transitions up to n = 31 have been calculated based on the 
simulation of a cold dark matter model and are given by Meiksin(2006). Inoue & Iwata (2008) adopt 
a different technique for calculating the absorption due to the Lya forest. However, both methods 
agree excellently for the intergalactic transmission at wavelengths greater than the Lyman limit (see 
Figure 6 in Inoue & Iwata 2008). Therefore, the treatment as prescribed by Meiksin(2006) is assumed 
throughout this investigation. Both models are based on recent estimates for the properties of the 
Lya forest, which show substantial differences from Madau (1995). 



Photoelectric absorption 

The contribution to the optical depth by photoelectric absorption for a source at redshift z is split 
into two parts (Meiksin 2006): the contribution from systems optically thin (tl < 1) at the Lyman 
limit and the contributions from the LLSs. The contribution from optically thin systems is revised 
from Meiksin (2006) according to 



rf ^ = A{1 + zl) 



4.4 



1 



1 + zl)^ (l + z) 



(2.2) 



where zl — — ^ and = 912 A. The normalisation is taken to be A = 0.07553, which is 

based on simulations. There is a degree of uncertainty in the amount of the absorption, depending 
on the ionization state of the hydrogen in the IGM and on the cosmological parameters assumed. 
The contribution due to the optically thick, tl > 1, Lyman Limit Systems is given by 



^LLS 



fdz' f 

JzL Jl 



OtlOz' 



1 — cxp 



1 + zl 
1 + z' 



(2.3) 



where q'^^q^, is the number of absorbers along the line of sight per unit redshift interval per unit 
optical depth of the system. The spatial distribution of LLSs adopted in Meiksin (2006) is 
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(2.4) 



3 



where A^n 



0.25 and 7 = 1.5. The optical depth distribution is given by 

dN _« 



drr 



OC Tr 



(2.5) 



where /3 = 1.5, although there may be small deviations from a perfect power law (Meiksin 2009). 

Inoue & Iwata (2008) treat the photoelectric contribution to the optical depth differently, as shown 
in Figure 6 of their paper. At wavelengths shorter than the Lyman limit there is some disagreement 
from the model used in Meiksin (2006). Inoue & Iwata treat the DLAs and LLSs separately, but they 
show that the contribution of DLAs to the Lyman continuum is small. They assume the distribution 
for tl given by Eq. (2.5) for the LLSs using (3 = 1.3, based on a fit by eye to recent observational 
data. A common z distribution for all types of absorbers is assumed, but this distribution normalised 
to only include the optically thick LLSs is given by 
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400. Figure [T| taken 



where they adopt z\ = 1.2, 
from Inoue & Iwata (2008), shows a comparison between the z distributions for the LLSs adopted 
by Meiksin (2006), Inoue & Iwata (2008), Madau (1995) and some recent observational data. It 
can be seen there is good agreement between the observed data and the distributions from the two 
more recent models. If the Inoue & Iwata (2008) option is selected, the code uses only their LLS 
distribution while retaining the model of Meiksin (2006) for the optically thin IGM and the Lya 
forest. 
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Figure 1: The number of LLSs per unit red shift along an average line of sight as a function of the LLSs' 
redshift. The filled circles are the recent observed data by Peroux et al.(2005). The open circles are the 
model assumed by Inoue & Iwata (2008). The dashed line is the assumed distribution in Meiksin (2006). 
The dotted line is that assumed by Madau (1995). 



The intergalactic transmission function 

The total effective optical depth of the IGM, Tes, is the sum of the resonance line and photoelectric 
contributions. It is a function of the observed wavelength, as well as the redshift of the source, and 
is given by 

reff = rLc + ri«^ + $^T-, (2.7) 
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where Tlc is the total Lyman continuum contribution from all LLSs, is the contribution from 

the optically thin IGM, Eq. (2.2), and is the optical depth for the ith Lyman series line. The 



attenuation of a background source by the IGM may then be characterised by the transmission 
function exp(— Teff). An example of the mean transmission through the IGM, taken from Meiksin 
(2006), is given in Figure [2] A value of 1 represents complete transmission through the IGM, at 
an observed wavelength Xohs, while a value of represents no transmission (hence the source is 
blacked out at these wavelengths). Each step observed in the transmission function, in the region of 
wavelengths greater than the Lyman Limit (Aobs > (1 + 2)912 A), represents the addition of a new 
line in the Lyman series. The Lyman continuum contributes at wavelengths shorter than the Lyman 
limit. 
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Figure 2: Intergalactic transmission as a function of observed wavelength. Mean intergalactic transmission 
for a source at redshift z = 3.5 (solid lines). The mean intergalactic transmission assuming no LLSs lie 
along the line of sight is shown by the dashed lines. (N.B. This uses an older model for the contribution 
to the optically thin IGM.) The thick lines are those using the model assumed by Meiksin(2006); the light 
lines are the estimates of Madau(1995). 



2.2 Broad band magnitudes, colours and IGM k-corrections 

The apparent AB-magnitude of a source of intrinsic flux f\ (in CGS units) measured through a filter 
with transmission T(A) and effective optical depth of attenuation. Teg, (see previous section) is given 

by 

mAB = -2.5\og,oj rflogA/A(AVc)e-^^«r(A) -48.59 (2.8) 

(Fukugita et al. 1996), where A is the wavelength. The code uses the normalisation convention for 
T(A) of 

POO 

/ dlogAT(A) = l. (2.9) 
Jo 



Alternative magnitudes on the Vega-magnitude system, for which the star Vega is defined as 
having zero magnitude in all bands, are also allowed for. 

The difference between the magnitudes with and without intergalactic attenuation is defined as 
the IGM /;;-correction 

kiGM = rriABircs) - mABir^s = 0), (2.10) 

where mAB(Tefr) is the magnitude through the IGM with effective optical depth t^s, and mAB(Teff = 0) 
is the magnitude if there were no attenuation due to the IGM. 
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3 Code features 



3.1 Calculating intergalactic transmission functions 

The code is designed to reproduce the intergalactic transmission curves in Meiksin (2006) by using 
Monte Carlo simulations. The average number N of LLSs is calculated by integrating the assumed 
distribution, Eq. (2.4), between z = and z (the chosen redshift of the source). The number 



of LLSs actually intercepted along each line of sight is then drawn from the Poisson distribution, 
f{k, N) = N^e~^ /k\, where k is the number of systems. For k > 10, the limiting Gaussian distribu- 
tion is used. For each LLS, a value of z and tl is assigned drawn randomly from the distributions Eq. 



(2.4) and Eq. (2.5). The contribution to attenuation by the diffuse IGM, Eq.(2.2), is also included. 



The normalisation A is an input parameter in the GUI. 



For each LLS, with optical depth and redshift z, 
depth at a given wavelength A is calculated using 



the photoelectric contribution to the optical 



tlc(A) = -Tl 



X 



Al(1 + ^) 



(3.1) 



which is derived from Eq. (2.3) and Xl = 912 A. This is calculated over a certain wavelength range 



with a given resolution. The default resolution is 1 A, but this may be altered by the user. The 
contributions from each LLS along the line of sight are then summed. The mean contribution due to 
each line in the Lyman series as derived in Meiksin (2006) are summed on top of the photoelectric 
contribution. 



The process is repeated for each realisation (i.e., line of sight) and an average is then taken to 
get the total effective optical depth of the IGM and hence transmission function exp(— Tofj) (see 
Section 2.1 ). 



An option is included to distribute the redshifts of the LLSs according to that of Inoue & Iwata 



(2008), Eq. (Q. As discussed in Section |2l| Meiksin (2006) and Inoue & Iwata (2008) agree closely 
on the resonance line scattering contribution to the IGM transmission. The user may use this option 
to examine the effect of distributing the LLSs according to the empirically derived z distribution 
of Inoue & Iwata, although the model used by Meiksin (2006) for resonance line scattering and the 



optically thin IGM (Eq. (2.2)) is retained. This option is implemented using the same method as for 
the Meiksin (2006) distribution. The parameters of the distribution function are included as input 
options for generality; these may be adjusted by clicking on the appropriate tab, which enables the 
user to enter alternative values. 

It is also possible to compute the effect of a fixed set of LLSs on the galaxy colours instead of 
using a Monte Carlo realisation. This may be done by adding a file in the top-level directory (such 
as input.dat) containing in the first line the number of LLSs, followed by a line for each absorber 
containing the data pair z and ti (with a space between), corresponding to the redshift and Lyman 
limit optical depth of the absorber. 

3.2 Calculating broad-band magnitudes, colours and IGM k-corrections 

The code computes the broad-band magnitudes and colours of model galaxy spectra for a choice of 
broad-band filters. 
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Filters 

A large choice of filter transmission curves (see Section 2.2) were included and can be selected from 
drop-down menus in the GUlj^ The following filters have been implemented: 

• Hubble Space Telescope (HST^ the filter responses for F300W, F450W and F606W for the 
Wide Field Planetary Camera 2 (WFPC2), referred to as f/soo, -B450 and Veoe; the filter re- 
sponses for F435W, F606W, F775W and F850LP for the Advanced Camera for Surveys (ACS), 
referred to as ^435 (approximating Johnson B), Vqqq, ijj^ (approximating SDSS i) and Z^^q (ap- 
proximating SDSS z)] the filter responses for FllOW and F160W for the Near Infrared Camera 
and Multi-Object Spectrometer (NICMOS), referred to as Jno and i^ieo- 

• UK Infrared Telescope (UKIRT): the filter responses for the UKIRT Infrared Deep Sky Survey 
(UKIDSS) filters if, J, K, Y and Z; the filter responses for the Wide Field Camera (WFCAM) 
filters J and K. 

• Subaru XMM-Newton Survey (SXDS^^ the filter responses for i?, V , R, i' and z' . 

• Keck: the filter responses for the Low Resolution Imaging Spectrometer (LRIS) B, V, Rg, Rkc 
(Kron-Cousins) and V. 

• Palomar 200-inch: the filter responses for f/„, G and R. 

• Canada-France-Hawaii Telescope ( CFHT): The filter response for u. 

The response curves also include estimates for the total system throughput (i.e., including the 



response of the full system). Filter transmission curves T(A) are normalized according to Eq. (2.9). 
Gaussian quadrature (see Appendix [B]) is used to compute the normalisation constant. 

The user should verify each filter response against expectation before performing any quantitative 
analysis. Additional filter response curves may be readily incorporated by adding files containing 
them to the filters directory with the following format: the first line should contain len and n 
(separated by a space), followed by len pairs of A and T(A) (separated by a space), where len is the 
number of wavelengths, n is a multiplyer of the wavelength A such that nX is the wavelength in A, 
and r(A) is the filter transmission at wavelength A. 

Spectra 

A variety of different intrinsic galaxy spectra are available for selection, taken from the STARBURST99 
model of Leitherer et al. (1999). The GUI is designed so that an age and metallicity of the model 
galaxy may be chosen, as well as an option to chose whether to use an instantaneous or continuous 
star formation model. The models use a Salpeter Initial Mass Functioij^ with stellar masses in the 
range 1 < M/Mq < 100. The models include emission from both stars and nebulae. The spectra 
may be chosen from a drop down menu and radio buttons. The fx values are linearly interpolated to 
the same wavelength values as the normalized transmission functions T(A) of the filters. The choice 
of redshift z for the model galaxy is included as an input parameter. The wavelengths of the spectra 
are consequently shifted by (1 + z). Internal reddening by gas and dust within the galaxies is not 
included. 



^This catalogue can easily be extended by the user, if a particular survey is to be investigated, as explained in 



Section 3.2 

^http://www. stsci.edu/hst 

lOprovided by M. Cirasuolo (private communication). 

^"'^The Salpeter Initial Mass Function is given by £,{M) = ^oM~'^-^^, where S^{M)dM is the number of stars born 
with masses between M and M + dM and ^0 is a normalisation constant. 
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Additional spectra may be readily incorporated by adding files containing them to the uploadedSpectra 
directory with the following format: the first line should contain len and n (separated by a space), fol- 
lowed by len pairs of A and fx (separated by a space), where len is the number of wavelengths, n is a 
multiplyer of the wavelength A such that nX is the wavelength in A, and fx is the wavelength-specific 
flux at wavelength A. 

Calculating magnitudes and colours 



The transmission curves exp(— Tgff) are calculated using the methods outlined in Section 3.1 and 
then linearly interpolated to the same wavelength values as the filter transmission curves T(A). 
The interpolated values of exp(— Tcfr), T(A) and redshifted fx are then multiplied together and the 



AB-magnitudes calculated using the integral (2.8), again with Gaussian quadrature. An option is 



provided to compute magnitudes referenced to Vega. 



IGM k-corrections 



The IGM k-correction, Eq. (2.10), for each filter is computed by also calculating the AB-magnitude 
with no intergalactic attenuation (i.e., Tes = 0) and subtracting this from the magnitude including 
attenuation. 

3.3 Distribution of absorbers 

Figure [3] shows the resultant z distribution of LLSs as produced by the Monte Carlo simulation. 
This test was done by running the simulation for a galaxy at redshift, z = 8, for 2000 realisations 
and outputting the properties of each LLS. These were put into redshift bins of width, Az, and then 
plotted using the approximation dN/dz ~ AA^/Az (averaged over all realisations), where AA^ is 
the number of absorbers produced in each bin. The figure shows the distribution of LLSs according 



to both Meiksin (2006) (Eq. and Inoue & Iwata (2008) (Eq. (Q), using the parameters 

suggested in their respective papers. The results of the Monte Carlo simulations are seen to match 
the expected distributions. The data are fairly noisy at the lower redshifts due to the small number 
of LLSs in each redshift bin. 

3.4 Calculating the mean intergalactic transmission functions 

Example mean transmission curves as produced by the Monte Carlo simulation are shown in Figure 
111 These are calculated assuming the model for the Lya forest as given in Meiksin (2006) and 



the contribution due to the optically thin IGM, according to Eq. (2.2). The result of distributing 



the LLSs according to both Eq. (2.4) and Eq. (2.6) are shown for comparison. The steps above 



the Lyman limit due to the Lyman transitions are observed. The Lyman continuum below the 
Lyman limit, due to the LLSs and optically thin IGM, can clearly be seen. The transmission in 
the continuum regime according to Inoue & Iwata (2008) is smaller because of the generally larger 
number of LLSs along each line of sight (as demonstrated in Figure [s]). 

3.5 Transmission along individual lines of sight 

There is a stochastic nature of the amount of attenuation along different lines of sight due to the 
rare LLSs. The result is that the transmission in the Lyman continuum regime (Aobs < (1 + 2^)912 
A) can vary dramatically along different lines of sight. A mean contribution from the Lya forest is 



assumed (see Section 2.1). Figure |5] demonstrates the transmission for three separate lines of sight 
for a source at z = 3. In the two extreme spectrum could be completely truncated below 

the Lyman limit, or if there are no LLSs near the source, significant transmission could still occur. 
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Figure 3: Number of LLSs per unit redshift along an average line of sight. The black points represent 
the output of the Monte Carlo simulation using the Meiksin(2006) model with 2000 realisations. The black 



dashed line shows the expected result using Eq. (2.4). The blue diamonds represent the output of the 



Monte Carlo simulation using the Inoue & Iwata (2008) model. The blue dashed line shows the expected 



results using Eq. (2.6) 
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Figure 4: Mean intergalactic transmission, over 2000 lines of sight, as a function of observed wavelength, 
Aobs) for sources at z = 3, 4 and 5 as viewed from left-to-right. The heavy lines assume the distribution of 



LLSs given in Eq. (2.4), the light lines assume that given in Eq. (2.6). The values of /3 in Eq. (2.5) used 



are 1.5 and 1.3, respectively. 



Indeed, for a source at redshift z = 3, assuming Eq. (2.4), the average number of intervening LLSs 



expected is 3.1. Therefore, from the Poisson distribution, the probability that no LLSs occur at all 
is = 4.5%. 
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Figure 5: Example transmission curves for a source at z = 3 along three separate lines of sight, assuming 



the z distribution of LLSs Eq. (2.4) and (3 = 1.5. Panel (a) shows the transmission when there are no LLSs 
near the source; (c) when there are a several absorbers near the source and (b) is an intermediate case. 

4 Acknowledging use of the code 

The use of IGMtransmission should be acknowledged using a statement like 

This paper computed IGM transmission values using IGMtransmissior|^ (Harrison, 
Meiksin & Stock 2011), based on the transmission curves of Meiksin (2006). 



"Available for download from 



http://code.google.eom/p/igintransmission 

If the LLS distribution from Inoue &: Iwata (2008) is used, a reference to their work should be 
added as well. Any work using the model galaxy spectra should reference Leitherer et al. (1999). 
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A The Graphical User Interface 
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Figure 6: The GUI produced to allow easy input of the numerous options and parameters associated 
with the program. 

Filter transmission functions may be chosen from the drop down menus. The program calculates 
two colours simultaneously; for example, if the program was in the set-up shown in Figure |6] it will 
calculate colours U300 — B^^q and 1^450 — Vgoe for the Hubble Space Telescope, Wide Field Planetary 
Camera 2 (WFPC2) including estimates for the total system throughput. These are taken from a 
library of transmission curves in a directory associated with the GUI. This catalogue is easily ex- 
tended by simply adding transmission curve files to the directory, as described in Section |3.2[ The 
choice of AB-magnitude system or Vega-system is specified by selecting the appropriate radio button. 
If a starburst model is to be used, chcking on the Starburst99 radio button will allow the age and 
metallicity of the model galaxy spectra to be chosen from drop down menus. Either an instantaneous 
or continuous star formation model may be chosen by selecting the appropriate radio button. The 
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relevant spectra is then selected from the catalogue. The redshift to place the model galaxy at is 
input into the text field. It is now possible to calculate the colour of the model galaxy, without any 
absorption, by simply hitting Calculate Colour. The colour will be printed to the chosen output file. 
Alternatively, the user may upload their own galaxy spectra by clicking on the appropriate radio 



button. The format for the spectra is described in Section 3.2 



To investigate the effects of intergalactic absorption, the user may choose whether to use the redshift 



distribution of Meiksin (2006) (Eq. (Q) or Inoue & Iwata (2008) (Eq. (Q) for the Lyman Limit 
Systems (LLSs) by selecting the appropriate radio button. The parameters for these distributions 
may be changed by clicking on the appropriate tab and entering the desired values into the text 
fields. The opacity values are taken from (Eq. ( 2.5[ )), with variable parameter (3. If either of these 
distributions are selected, the program will perform Monte Carlo simulations for the number of lines- 
of-sight (realisations) defined, by randomly drawing opacity and redshift values for the LLSs from 
the chosen distribution. Alternatively, there is an option to use custom defined properties of Lyman 
Limit Systems from an input file, as described in Section [STTj If any of these methods are chosen, the 



contribution to the absorption by the Lya forest will be included according to the model of Meiksin 



(2006) and the contribution by the optically thin IGM according to Eq. (2.2) for which the normali- 



sation constant A may be chosen. The resolution for calculating the transmission may also be chosen. 



Several options are available for the output. By checking or un-checking the relevant box, it is 
possible to either print the colours for each realisation separately or print the average over all lines of 
sight to the user-defined output file. There is also an option to print to file the cumulative distribu- 
tion of colours. If absorption is included, the program will also print to file; the average transmission 
curve (see Section 3.4), the number and properties of the LLSs used for each line of sight and the 



average IGM k-correction values for each filter (as calculated from Eq. (2.10)), which may be useful 
for interpreting the results. 



B Integration using Gaussian quadrature 

When numerical integration is required, Gaussian quadrature has been used. This involves approx- 
imating the integral of a function f{x) by 

N 



I W{x)f{x)dx = J2^,f{x,) (B.l) 



where W{x) is the weight function, wj are the weights and xj are abscissa points. The function 
is evaluated at each xj value. The code uses the weight function W{x) = 1, corresponding to the 
Gauss-Legendre case. The weights in this case are given by 

2 

where are the Legendre Polynomials. A routine to compute the integral between xi and X2 was 
adapted from Numerical Recipes in C (Press et al. 1992), and translated into Java. The transmission 
functions T{X), the galaxy spectra f\ and transmission functions exp(— r) were all evaluated at the 
same absicissa points using linear interpolation. This method was tested using a simple Simpson's 
rule routine to verify that both methods agreed. 
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